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Abstract 

In this article we introduce a new method for constructing implicit 
symplectic maps using special symplectic manifolds and Liouvillian forms. 
This method extends, in a natural way, the method of generating functions 
to 1-forms which are globally defined on the symplectic manifold. The 
maps constructed by this method, are related to the symplectic Cayley’s 
transformation and belong to a continuous space of dimension n(2n+l). 
Applying the implicit map to the discrete Hamilton equations we obtain 
the generalized symplectic Euler scheme. We show the relations of the 
elements of this family with other discrete symplectic mapping, in par¬ 
ticular 1) with the mappings obtained by generating functions of type 
I, II, and III and IV; 2) with the symplectic Euler methods A and B; 
and 3) with the mid-point rule. Moreover, we show the corresponding 
symplectic diffeomorphisms and their Liouvillian forms on the product 
symplectic manifold. We illustrate the details of the method in construct¬ 
ing two different families of implicit symplectic maps for n = 1. This is 
a geometrical method which overcomes the difficulties of the Hamilton- 
Jacobi theory and generating functions. 


1 Introduction 

Symplectic maps can be constructed by using generating functions, and they 
were introduced by Poincare when looking for periodic orbits of second genus 
m- From the numerical point of view, symplectic maps are used for simulating 
Hamiltonian dynamics, based on the fact that the Hamiltonian flow is a one 
parameter subgroup of symplectic diffeomorphisms. In the classical construc¬ 
tion of symplectic maps, the use of Darboux’s coordinates in different stages 
of the procedure hides several interesting properties, which arise with the dis¬ 
cretization in time of Hamiltonian flows. Other properties do not appear in the 
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continuous case due to a classical procedure of reduction from symplectic to 
contact geometry, which takes the Hamiltonian vector field into the Reeb field 
on the reduced manifold 00- 

In this paper, we address the problem of constructing implicit symplectic 
maps with a different approach. Instead of using Hamilton-Jacobi theory to 
obtain a generating function, we take advantage of the geometrical properties 
of both the Hamiltonian flow and the product symplectic manifold. However, 
we avoid the use of Darboux’s theorem and the classical evolutive Hamilton- 
Jacobi equation. The reason is that the coordinates of the original problem are 
locally conjugated to the canonical coordinates by a symplectomorphism which 
is hidden by Darboux’s theorem. In consequence, the classical construction 
integrates a modified problem in normal form using action-angle coordinates. 
Instead of Darboux’s theorem we use special symplectic manifolds mmm 
and their Liouvillian forms in order to track the state of the original coordinates 
at every stage of the procedure. 

This approach differs from the standard procedure of generating functions 
in two ways. First, we consider Lagrangian submanifolds as the integral sub¬ 
manifolds of a prescribed Liouvillian form 9 00, i.e. submanifolds defined by 
the differential ideal I generated by 9 0. Second, we use the complex struc¬ 
ture on the product manifold to obtain the equations of the tangent space to 
the Lagrangian submanifold defined implicitly by the Liouvillian form. The 
advantage in the symplectic case is that the vertical bundle corresponds to the 
characteristic bundle to the Lagrangian submanifold. In our method, the Liou¬ 
villian form substitutes the generating function and the multiplication with the 
complex structure on the product manifold substitutes the role of the resolution 
of the Hamilton-Jacobi equation. An additional projection of the Lagrangian 
submanifold offers the possibility to construct symplectic maps by a simple 
transversality condition. For the linear case, this condition is weakened by the 
classical condition stating that the symplectic map must converge to the identity 
map on the diagonal of the product manifold. All this procedure is constructed 
on the framework of special symplectic manifolds. 

The latter condition is translated into a condition relating two different 
Liouvillian forms on the original symplectic manifold. Surprisingly, using this 
relation, the framework of special symplectic manifolds is no more necessary. In 
this way, we can construct symplectic maps just using Liouvillian forms, leading 
to the method of Liouvillian forms. This procedure reveals some remarkable 
properties and gives a different perspective on Liouvillian forms and generating 
functions. We give an interpretation of this method as a discrete version of 
the Cayley’s transformation of a Hamiltonian matrix, leading naturally to a 
symplectic map. Finally, we construct a numerical scheme, which is the natural 
generalization of the Euler symplectic schemes. 

We summarize the rest of the paper. In Section 2 we give the main definitions 
and classical results on symplectic manifolds and the product manifold that 
are necessary for our exposition. Special symplectic manifolds and Liouvillian 
forms are exposed in Section 3. In this section we introduce the generating 
functions of type I, II, III and IV and their related Liouvillian forms expressed 
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as elements of some special symplectic manifolds. We return systematically to 
these examples in the rest of the paper. In Section 4 we expose the way we can 
construct symplectic maps using Liouvillian forms. Section 5 gives the structure 
of Liouvillian forms on the symplectic manifold and the product symplectic 
manifold. In Section 6 we construct the implicit map. Section 7 describes all 
four steps of the method of Liouvillian forms and we state our main result. In 
addition, we relate these maps to the symplectic Cayley’s transformation. In 
Section 8 we apply the implicit maps to Hamiltonian dynamics and we introduce 
the generalized implicit symplectic Euler scheme. We show that the mid-point 
rule, and the symplectic Euler A and B maps are particular cases of our scheme. 
Finally, in Section 9 we develop two examples of continuous families of implicit 
symplectic integrators: the first one using special symplectic manifolds and the 
second one applying the generalized symplectic Euler scheme. 


2 Symplectic manifolds and symplectic mappings 

In this section we recall some classical results and definitions in order to uni- 
formize the notation. The results are stated in a geometrical fashion in prepa¬ 
ration of the next section where we will set the framework of the method of 
Liouvillian forms. We assume the reader is familiar with the terminology of dif¬ 
ferential geometry and vector bundles. For an introduction the reader is referred 

to mm i- 

A symplectic manifold is a 2n-dimensional manifold M equipped with a 
non-degenerated, skew-symmetric, closed 2-form w, such that at every point 
m £ M, the tangent space T m M has the structure of a symplectic vector space. 
In addition, we say that (M, to) is an exact symplectic manifold if there exists 
a 1-form 9 such that oj = d9 ; 6 is called a Liouvillian form. In what follows, all 
the symplectic manifolds will be exact. 

An almost complex stucture J on a manifold M is a smooth tangent bundle 
isomorphism J : TM —>• TM covering the identity map on M such that for each 
point z £ M, the map Jz = J{z) : T Z M —> T Z M is a complex structure on 
the vector space T Z M , it means an endomorphism on T Z M such that J' 2 (v) = 

J o J{v) = — v for every v £ T Z M. We write J 2 = — I for simplicity. A 
symplectic manifold with an almost complex structure possesses a Riemannian 
structure g which enables the definition of a symmetric positive definite form 
where <*;(•,•) = <?(•, J-). We fix the Riemannian structure = (•, •) such 

that at every point on any symplectic manifold x £ M, the tangent space 
( T X M , (•, -) x ) is an inner product space. 

A submanifold A C M is called Lagrangian if the symplectic form restricted 
to A vanish w|a = 0. The tangent bundle to the symplectic manifold TM splits 
into the tangent T A and the vertical VA bundles of the Lagrangian submani¬ 
fold in the form TM = T A ® V A, where VA := (TA) X is the complementary 
subbundle to TA. The complex structure J and the non degeneracy of to on M, 
let us rewrite the vertical bundle as V X A = JT X A for every x £ A. We write 
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this as 


T X M = T X A ® J{T X A). 

Let (M, w) be a symplectic manifold. A diffeomorphism / : M —> M is said 
symplectic or a symplectomorphism if f*u) = ui, where /* denotes the pullback 
of the 2 -form defined by 

(f*u) z {u,v) = w f ( z )(Tf(u),Tf(v)), z,f(z) G M,u,v eT.M. (1) 

In this expression Tf:TM —> TM is the tangent map of / : M —> M. The set 
of all the symplectic diffeomorphisms on a symplectic manifold form a group 
denoted by Sp(M,w) and called the symplectic group of (M,w). 

In order to construct symplectic maps we follow the classical construction. 
Define the product manifold of two copies of an exact symplectic manifold (M, ui) 
at times t = 0 and t = h, which we denote by (Mi,wi) and (M 2 ,w 2 ), respec¬ 
tively. The manifold M = Mi x M 2 with the canonical projections 7 r,; : M —> Mj 
for i = 1.2 let us define the forms Oq and ojq on M by 

9q = 7T*6»i - 71-202, W 0 = ITiUl! - 7rjw 2 (2) 

We have the following facts (see [U sec 5.2] for the proofs): 

• (M,ujq) is a symplectic manifold of dimension An. 

• For any symplectic map 4> '■ Afi —>• M 2 , the graph of denoted by 
and defined by 

A^ = | (Z) e M | z e Mi,(f>(z) G M 2 | 

is a Lagrangian submanifold of M. It means Wq|a^ = 0 

• Since ujq = cIQq, then 9q is a locally closed form on A^. Applying 
Poincare’s lemma, 6q is locally exact in a neighborhood of A^. Conse¬ 
quently, there exists a function S defined on the Lagrangian submanifold 
A^ such that its differential concides with the restriction of the 1 -form 9q 
to A 0 . It means that 

dS\\+ — |a* (3) 

The function S : A^ —> R is called a generating function for the symplectic 
map </>; 

• There exists and induced endomorphism on T X AI which becomes the as¬ 
sociated complex structure on w e given by 

J = (Tmf Ji - (Ttt 2 ) t J 2 


where Ji, are the associated complex structures to u>i, i = 1,2. 
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Symplectic maps, generating functions, Liouvillian forms and Lagrangian 
submanifolds are closely related. For instance, in a generic symplectic manifold 
M, any Lagrangian submanifold Ac M which is transverse to the fibers of the 
projection ttm '■ T*M —> M can be parameterized by a suitable (local) function 
S. In contrast, Liouvillian forms are globally defined and they do not need a 
particular parameterization, however topological properties of the underlying 
manifold like characteristic classes prevent the use of these forms in particular 
applications. 

We rephrase some classical properties of Lagrangian and symplectic sub¬ 
manifolds in term of the symplectic product manifold M. 

Lemma 2.1 Let A C M be a Lagrangian submanifold and $ £ Sp{M, u>q) a 
symplectomorphism. We have the following facts: 

1. The image of the Lagrangian submanifold under $ is again a Lagrangian 
submanifold of M. 

2. The projection 7Tj(A) C Mj is a Lagrangian submanifold in Mi, i = 1,2. 

Proof. 1) Let A = <h(A) be the image of A under the symplectomorphism <f>. For 
all A and £, 77 G T y A there exist x G A and u, v G T x A such that y = $(x), 
£ = T$(u) and rj = T$(i>). Since = lo q we have 

(ve) y (€,v) = (u e ) x (u,v) = 0 

and A is a Lagrangian submanifold in (M,cuq). 

2) Since M = M\ x M 2 is a product manifold, its tangent bundle splits into 
TM = TM\ ® TM 2 . The subblundle TA C TM is given under the splitting 
by TA = TAi © TA 2 . We have T7r,;(TA) = TA,; c TM.i = 1,2. Due to 
the spitting TA = TAi © TA 2 then Wq|a = 0 if and only if = 0 and 

consequently 7r,;(A) C A/,; are Lagrangian submanifolds for i = 1, 2. □ 

Lemma 2.2 Let <f>i G Sp(Mi,u>i), i = 1,2 be two symplectomorphisms. The 
induced diffeomorphism on M by diagonal action on M±xM 2 defined by <F(M) = 
4>i(Mi) x (j> 2 (M 2 ), is a symplectic diffeomorphism in M. 

Proof. It is enough to show that >F*u;e =Uq. By definition of the symplectic 
form and the fact that <f>i, i = 1 , 2 are symplectomorphisms we have 

UJq = 7T*CCi — 7^2002 = O (fr^COi — 7r 2 ° 0 2 CU 2 = (01 O 7Tx)*CCi — (<(> 2 O 7T 2 )*CU 2 . 

Applying the properties of the pull-back of the composition we obtain 

COq = (01 O 7Tl)*Wl — (0 2 O 7T 2 )*W 2 = $*W e 

which gives the result. □ 

This last result gives us the possibility to consider Lagrangian submanifolds 
which do not look like the graph of a diffeomorphism in M. It lets us work with 
Lagrangian submanifolds in mixed coordinates, or in other words with generic 
implicit symplectic mappings. 
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3 Special symplectic manifolds and Liouvillian 
forms 

An important class of symplectic manifolds are the cotangent bundle to Rieman- 
nian manifolds. In particular, they model the phase space of many Hamiltonian 
mechanical systems, where the base coordinates correspond to positions or con¬ 
figurations, and vertical coordinates correspond to momenta. 

Let Q be a C°° manifold. Consider the cotangent and tangent bundles 
T*Q and TQ, and its canonical projections on the base manifold Q denoted by 
tq-.TQ^Q and -kq : T* Q —»• Q respectively. 

The projection ttq determines a natural map between the cotangent bundle 
and the double cotangent bundle by its pullback ttq : T*Q — >• T*(T*Q), which 
sends 1-forms on Q to 1-forms on T*Q by 9 = ir*r] for every 1-form ?y defined 
on Q. Formally p is a section of the cotangent bundle that we denote by 77 S 
T(T*QL w j iere y(T*Q) denotes the space of smooth sections on the cotangent 
bundle @ 


9 Gr T*(T*Q) 



The form 9 = 9 V , induced by the identity morphism 9 v ( 77 ) for every 

rj G T(T*Q), is called the Liouville form on T*Q and is alternatively defined by 
its action on the tangent bundle by the equation 

(v,9 v ) t *q = (Tn Q (v),r 1 ) Q , v G T V T*Q, r/ G T(T*Q), 9 V G T*(T*Q). 


v G TT*Q -(vine- T*(T*Q) 3 9„ 



TQ -- T*Q 9 r V 



Remark 1 The double cotangent bundle (T*(T*Q),T*Q,ttt*q) is defined by 
the projection 7 tt*q '■ T*(T*Q) —> T*Q. Consequently, the canonical Liouville 

1 We use the notation (j Gr E instead of /3 G r( E) for simplifying the notation in the 
diagrams. 
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form 9 £ r (T*(T*Q)) is a section in T*(T*Q) corresponding to the inclusion 
of sections from F(T*Q) into T(T*(T*Q)). This inclusion is also interpreted 
as the identity map. The fact that the Liouville form is a section of the double 
cotangent bundle is not evident when we work with symplectic vector spaces 
and it has been the source of many misunderstandings in some areas like the 
numerical analysis. 

When we shall have the occasion to deal with cotangent bundles of different 
manifolds, we will denote the Liouville form on T*Af by Ojy. 

The cotangent bundle T*Q inherits a natural symplectic structure to = d0Q 
such that the couple (T*Q,w) becomes a symplectic manifold. Unfortunately, 
many geometrical properties of symplectic manifolds lost significance in me¬ 
chanical systems due to a missing physical interpretation. Tulczyjew proposed 
in the study of symplectic manifolds which are diffeomorphic to cotangent 
bundles by means of special symplectic manifolds (see also mm)- 

A special symplectic manifold is a quintuple (M, Q, 9 , 7 r, ip) where n : M Q 
is a fibre bundle, 9 is a 1-form on M and p : M —> T*Q is a symplectic 
diffeomorphism such that 7r = 7rg o and 9 = <p*9q. 

9 £ r T*M ^ -—- T*(T*Q) 3 r 9 Q 



The pair ( M,d9 ) is a symplectic manifold. We call the 1-form 9 = <p*9q 
a Liouvillian for on M. More generally we say that a 1-form 77 on an even 
dimensional symplectic manifold (M, ui) is Liouvillian or of Liouvillian type if dr/ 
is a symplectic form on M, or equivalently if (M, dij) is a symplectic manifold. 

Let K C Q be a submanifold and S : K — y R a function on K. The 
Lagrangian submanifold generated by S in the manifold (M, d9) is defined by 
the equation (v,9) = (Tn(v),dS) where v £ TM and tm(v) = to, in the 
following way 


A = {m £ M|7 t(to) £ K, {v , 9) = ( Tn(v),dS }} . (4) 

We suppose K C Q is an embedded submanifold and we write it as an 
inclusion * : K Q. In this case i*(dS) £ T(T*K) and dS £ T(T*Q). If in 
addition, we suppose that i : K Q is an immersion, 7r defines a submersion 

2 Libermann and Marie [8] and Libermann [?j define Liouvillian forms in the more gen¬ 
eral context of bundle morphisms. We use a simplified definition for working with special 
symplectic manifolds. 
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and the submanifold A is well defined. The following graph shows this situation 


v G TM T*M 9 r 9 



R 


Any special symplectic manifold (M, Q, 9, 7r, determines a natural split¬ 
ting of the tangent bundle TM = T A© VA using the Liouvillian form 9 = ip*9 q 
and the projection 7r = ttq o <p. The projection keeps track of the deforma¬ 
tion of the tangent (horizontal) bundle, and the Liouvillian form recovers the 
deformation of the vertical bundle. 

The following result relates the symplectic manifolds of interest in this work: 
(M,w e ) and (T*(Qi x Q 2 ), d6 QlX Q 2 ), where 9q 1 x q 2 is the Liouville form on 
Qi x Q 2 • 

Proposition 1 Define coordinates (q,p,Q,P) G M in the product manifold 
such that ( q,p ) G Mi and ( Q,P ) € M 2 - Let E\ : M —>■ T*(Qi x Q 2 ) be the 
linear map given by 

Ei : M -»• T*(Qr x Q 2 ) 

( q,P,Q,P ) ( q,-Q,p,P ). (5) 

ITien is a symplectomorphism. 

Proof. Define coordinates (x,X,y,Y) g T*(Qi x Q 2 ) where (x,X) £ Qi x Q 2 . 
We have (x, X, y, Y ) = (g, —Q,p, P) and the Liouville form in these coordinates 
becomes 9q iX q 2 = ydx + YdX. Taking the differential we have 

E* (dy A dx + dY A dX) = dp A dq — dP A dQ = w@ 

as we want to show. 

□ 

Ei is known as the canonical symplectomorphism between the product man¬ 
ifold (M,w q) and the contangent bundle (T*(Qi x Q 2 ),d9Q lX Q 2 ). 

Corollary 1 Let G Sp(M ,Wq) be a symplectomorphism on M and consider 

Ei as below. Then the diffeomorphism T : M —> T*(Qi x Q 2 ) given by 


T = Ei o $ 


(6) 



is symplectic. 

Proof. Direct from the fact that ft : M —> T*(Qi x Q 2 ) is the composition of 
two symplectomorphisms 


M -^-«- M 



T*(Qi x Q 2 ). 


□ 

The quintuple (M, Qi x Q 2 ,9,n,'i t), where ft is defined in (JG]) becomes a 
special symplectic manifold on Qi x Q. 2 - 


T*M 


M 


7 r QlxQ 2 

Qi x Q 2 

We recall that any Lagrangian submanifold A £ T*(Qi x Q 2 ) defines a 
Lagrangian submanifold in M by the symplectomorpliism ’t” 1 (A). 

3.1 Examples of special symplectic manifolds 

We consider some simple symplectomorphisms between the product of two 
symplectic manifolds and the cotanget to the product of two configuration 
spaces M —>• T*(Q 2 x Q-f). These symplectomorphisms define four special 
symplectic manifolds and we will be interested in their Liouvillian forms. Us¬ 
ing symplectic coordinates on the manifolds (xi, X i: yi,Yi) £ T*(Qi x Q 2 ) and 
(qi,Pi,Qi, Pi) £ M we consider four symplectic diffeomorphisms 

1. Ei(qi,pi,Qi, Pi) = (qi,—Qi,pi, Pi) corresponding to the canonical sym¬ 
plectic diffeomorphism from M —>• T*(Q 2 x Qi). The Liouville form on 
T*(Q 2 x Qi) corresponding to Oq iX q 2 is pulled back to the Liouvillian 
form 9 = ( E*)9q iX q 2 and it corresponds to the form 9 e on M. For 
instance, the Liouvillian form and the projection are 

0 e = pdqi - PdQi and ^{qi.Pi, Qi,Pi) = (qi, -Qi)- 

2. ft \i(qi,Pi,Qi,Pi) = ( qi,Pi,Pi,Qi)■ This diffeomorphism gives the Liouvil¬ 
lian form and projection given by 

9 = pidqi + QidPi and ^(qi,Pi,Qi,Pi) = (qi,Pi)- 




■ T*T*(Qi x Q 2 ) 


v T*(.Q 1 xQ 2 ) 


T*(Q 1 x Q 2 ) 
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3. ^ m {qi,Pi,QiiPi) = ( Qi,Pi,~Pi,~qi )• This diffeomorphism gives the fol¬ 
lowing Liouvillian form and projection 

9 = -qidpi - PidQi and ^{quPuQuPi) = (Qi,Pi)- 

4. 'i lv {qi,Pi,Qi,Pi ) = ( ~Pi,Pi,qi,Qi ). This diffeomorphism gives the fol¬ 
lowing Liouvillian form and projection 

9 = -qidpi + QidPi, and n(q tl p l: Qi, Pi) = (~pi,Pi). 

Since the Liouvillian form on M corresponds to a class of symplectomor- 
phisrns modulo symplectic rotations we need both elements, the Liouvillian form 
and the projection for fixing the symplectomorphism between these manifolds. 

The symplectomorphims E\, JL,, \I/j„ will be revisited in the next section 
where we related them to generating functions. 


4 Symplectic maps from Liouvillian forms 

The usual way to construct symplectic maps on the product manifold (M,Wq) 
is using generating functions S : A —> R defined on some Lagrangian sub¬ 
manifold A C M; it is an inverse problem. This inverse problem is solved using 
Hamilton-Jacobi theory for estimating the characteristic bundle which contains, 
as a subbundle, the vertical bundle to the Lagrangian submanifold. We deal 
with this problem in a more direct way using Liouvillian forms which define 
Lagrangian submanifolds as their integral submanifolds. The transformation 
between the vertical and the tangent bundles is given by the complex structure 
associated to the symplectic form. In this way we avoid the solution of both the 
Hamilton-Jacobi equation and the generating function. However, in both pro¬ 
cedures, we must select what type of symplectic maps we are looking for. In our 
case, we are interested on symplectic maps adapted for constructing numerical 
schemes which differ from those used for studying periodic orbits in an essen¬ 
tial way. The task in this section is to characterize Lagrangian submanifolds 
adapted for constructing numerical schemes. 

Remark 2 Symplectic maps for numerical schemes and those used for study¬ 
ing periodic orbits solve variational problems with different boundary conditions. 
The former consider the minimization of the action integral along paths join- 
ning two different fixed points; the latter consider closed paths with prescribed 
period T > 0. This implies that Poincare’s differential /onrd introduced in m 
for studying bifurcations of periodic orbits is not well suited for constructing 
numerical schemes. A detailed study of this fact is given in W- 

To obtain a symplectic map for constructing numerical schemes we need to 
recover a symplectic vector space from the information encoded in the tangent 

3 It is the differential of the so called Poincare’s generating function 
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space to the Lagrangian submanifold A C M. In fact, we state a stronger 
condition. Let (M, Qi x <2 2 ,#,7r, T) be the special symplectic manifold defined 
in the last section and consider the image N = n(M) of the projection 7r = 
7r CixQ 2 o^asa submanifold by the adapted inclusion i : N^-M. If we can give 
a symplectic structure u>n on N such that Wq = (7r*)w/v, and ojn = (**) w e then 
the integral submanifolds of the Liouvillian form 9 on (M,wq) are adapted for 
the construction of symplectic integrators. 

If the symplectic submanifold N C M belongs to a symplectic path joining 
z(t) £ Mi x Mi to z(t + h ) £ Mi x M 2 , then the submanifold (N,ojn) must 
converge to the original symplectic manifold 


lim (TV, un) = (Mi, wi) = (M 2 ,w 2 ). (7) 

h —^0 

We summarize these requirements in the following: 

Condition 1 Consider the tangent space to the projected manifold at a point 
n(x) £ N, x € M, given by T V ^N = Tn(T x M). Denote by x £ A C M a 
point on the diagonal of the product manifold. If 

1 . the restriction of w@ to the submanifold N = i r(M) is a symplectic form 
on N, and, 

2. Tn (t x M^ 

the Liouvillian form 6 of the special symplectic manifold (M, Qi x Q 2 ,(9,7r, 'L) 
defines a Lagrangian submanifold adapted for the construction of a symplectic 
scheme. 


= id, is the identity, 

.T —X 


Instead of recovering the symplectomorphism 'L which defines the special 
symplectic manifold, we use the fact that Lagrangian submanifolds are the 
integral submanifolds to Liouvillian forms. Let 6 i be a Liouvillian form on 
the 2n-dimensional symplectic manifold (M, w). For every symplectomorphism 
<j> : M —» M, the form 0 2 given by the pullback 6 i = < $*62 is again Liouvillian 
on (M, w). They are different Liouvillian forms producing the same symplectic 
structure u> = d 6 i = d0 2 on M. On the product manifold (M,uiq), the form 
9 = tt^Oi — it 2^2 is Liouvillian on M since dO = w e . 

Using point 2) in Lemma 12.11 and Condition [1] we will induce the Liou¬ 
villian form on the product manifold (M,uj q) by Liouvillian forms 9 1 and 0 2 
on the original symplectic manifold ( M , w) defining complementary Lagrangian 
submanifolds. For this, we consider the complex structure J associated to ui. 
There exists a (tautological) symplectomorphism J : M —>• M such that its 
tangent map is exactly the complex structure TJ = J : TM —> TM. This 
symplectomorphism is attached to every oj. 

Theorem 4.1 For every exact symplectic manifold ( M,ui), there exists at least 
a Liouvillian form 9 1 on the original symplectic manifold ( M , w) such that the 
Liovillian form 

0 = nlOi - n* 2 9 2 , 9 1 =r(9 2 ), (8) 
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on the product manifold ( M,uiq) satisfies Condition 1. 

We need additional elements for proving this theorem which follows from Lemma 
!5.2l below. For the moment, we will relate the Liouvillian forms of the examples 
given in section O with the classical generating functions. 

Lemma 4.2 The Liouvillian form of the symplectomorphisms E\, dbj, 4' m , 4b„ 
given in section \3.1\ are associated to the generating functions of type I, II, III 
and IV, respectively. 

Proof. We perform the same computations for every symplectomorpliism using 
the Liouvillian form 9 on the product manifold ( M, w e ) and the projection 
7r = tt(QixQ 2 ) ° ^ f° r the corresponding T. 

1. Ei produces the Liouvillian form 9 = pidqi — PidQi, which is locally equiv¬ 
alent to the differential of a function S : M —>• R with S = S(qi,Qi ) which 
is a generating function of type I. It Defines a Lagrangian submanifold in 
M by 

A = {{quPi, Qi, Pi) G M\Pi = §§r,Pi = —§&}■ 

It does not produces a map adapted for construct symplectic integrators 
since the projection Tr(qi,Pi, Qi, Pf) = ( qi,—Qi) does not give symplectic 
coordinates. 

2. dbj gives the Liouvillian form 9 = pidqi + QidPi, which is locally equivalent 
to the differential of a function S = S(qi,Pi ), which is of type II. It defines 
a Lagrangian submanifold in M by 

A** = {{qi,Pi,Qi,Pi) G M | pi = = ||-j . 

This Liouvillian form is adapted for constructing a symplectic integrator 
since the projection n(qi,pi, Qi, Pi) = ( qi,Pi) gives symplectic coordinates 
in the diagonal. 

3. Tjjj produces the Liouvillian form 9 = —qidpi — PidQi which is locally 
equivalent to the differential of a function S = S(Qi,Pi), which is of type 
III. It defines a Lagrangian submanifold in M by 

A ra = {(qi,Pi,Qi,Pi) G M\ qi = -§^,Pi = • 

This Liouvillian form is adapted for constructing a symplectic integrator 
since the projection tt( qi,Pi, Qi,Pf) = ( Qi,Pi ) gives symplectic coordinates 
in the diagonal. 

4. gives the Liouvillian form 9 = —qidpi + QidPt which is locally equiv¬ 
alent to the differential of a function S = S{pi, Pi), which is of type IV. It 
defines a Lagrangian submanifold in M by 

A«J = ^ {qi,Pi, Qi, Pi) G M\ qi = , Qi = ~ ~ 0 p[ ^ 4 
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This Liouvillian form is not adapted for constructing a synrplectic in¬ 
tegrator since the projection n(qi,pi,Qi, Pi) = (~Pi,Pi) does not give 
symplectic coordinates. 


5 Structure of Liouvillian forms 

Generating functions for constructing symplectic maps are those of type II and 
III, which is not the case for generating functions of type I and IV. This happens 
in the same way for Liouvillian forms. In this section we study the structure 
of Liouvillian forms for characterizing those adapted for constructing implicit 
symplectic integrators. We recall some facts about the Liouvillian forms on a 
generic exact symplectic manifold (M, uj). 


5.1 Liouvillian forms on an exact symplectic manifold 

Consider an exact symplectic manifold ( M,cL0 ), where 9 is a generic Liouvillian 
form. For every function F : M —> R. on M the differential form 9 + dF is 
again a Liouvillian form since d(9 + dF) = oj. It means, the set of Liouvillian 
forms is invariant under the addition of exact 1-forms. This comes from the fact 
that d?{F) corresponds to a symmetric tensor S = Hess(F 1 ). This generic fact 
implies that the space of Liouvillian forms is infinite-dimensional. 

For studying the structure of Liouvillian forms, we want to understand the 
structure of the symplectic form in a generic basis. Given local coordinates 
2 = { z i}i= l £ Af, a basis of T*M = u) is given by dz = {dzi}?^, and 

r \ 2 n 

a basis of T Z M = A 1 (M, a;) is given by d z = |^- j . These are dual bases 

satisfying dzi (^ 7 ") = ff 1 = &ij- I n th e basis {dzi}f n , the symplectic form to is 
given by 

2 n 

u> = y] \Jijdzj A dzi, Jij G M. (9) 


The matrix J = [Jij] satisfies J 2 = — I^n, and it is the representation of the 
complex structure associated to u in the basis {dzi}f n . Remark that J is an 
antisymmetric matrix J + J T = 02 n - 

A Liouvillian form 9 on ( M,u >) is expressed in the basis {d 2 i}?=i by 9 = 
cxi(z)dzi z = (zi, Z2, ■ ■ ■ , Z2n) G M, where on : M — > R are smooth functions 
of z. Since 




daj(z) 

dzi 


dZi A dzj, Z = (zi,Z2,---,Z2n) G M, 


and writing A(z) = [Aij{z)\ = 
and antisymmetric components by 


dai(z ) 
dzj 


we have the decomposition in symmetric 


A s (z) = \ (A(z) + A T (z)) and A a (z) = \ (A(z) - A T {z)) . 
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The condition for 9 to be a Liouvillian form is A a (z) = ij, and denoting 
A s (z) = S(z) to remark it is a symmetric matrix, we have A(z) = S(z) + \J. 
Moreover, the symmetric part S(z) belongs to the kernel of the differential 

9 A- to. We have proved the following: 

Lemma 5.1 The set of Liouvillian forms on an exact symplectic manifold (M, w) 
is given in local coordinates by 1-forms 9 = 'f2 i ai{z)dzi where the matrix 
A{z) = da Q^ has the structure A(z) = S(z) + \ J, where S(z) = S T (z) and J 
is the complex structure associated to oj in the local coordinates {zi}f n . 

Every exact symplectic manifold ( M , u>) possesses a natural Liouvillian form 
$o given by the complex structure associated to the symplectic form w. The 
Liouvillian form 9q has no symmetric component and we call it the elementary 
or basic Liovillian form. 

Corollary 2 The elementary Liouvillian form in Darboux’s coordinates z = 

( q,p ) £ (M,lo) is given by 


do = \ ^2 (Pi d 1i ~ 9idpi ) ■ (10) 

2 = 1 

Consider the tautological symplectomorphism J as in Theorem 14. II 

Lemma 5.2 The elementary Liouvillian form 9 q is invariant under the action 
of the complex structure, i.e. 9 q = J*9q. 

Proof. Given local coordinates z £ (A 1, oj) the elementary Liouvillian form is 
given by 9 o = \dzJz. Consequently 

3*9 0 = \{dzJ T )J(Jz) = \dzJz = 9 0 . (11) 

□ 

This lemma gives a Liouvillian form on an exact symplectic manifold, prov¬ 
ing Theorem 14.11 The Liouvillian form 9 = 7r*(?o — 9q induces a natural 
symplectomorphism which will be associated to the mid-point rule. We will 
show this fact in the last section where we will construct numerical integrators 
from the symplectic maps coming from Liouvillian forms. 

For studying the way a Liouvillian form induces a symplectic integrator we 
restrict our study to the set of Liouvillian forms with linear coefficients. It 
becomes a linear space over R with finite dimension. In the rest of this paper 
we replace the matrix A(z) for a linear form Az = (S + \ J)z where A and S 
are constant matrices in M 2 , 2 x 2 ra(R)- 

Lemma 5.3 The space of Liouvillian forms with linear coefficients on a 2n- 
dimensional exact symplectic manifold has dimension n(2n+ 1). 

4 Remark that dz is a covector and it transforms by the transpose dz H> dzJ T . 
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Proof. Applying Lemma 15.11 we deduce that the dimension of the space of 
Liouvillian forms is exactly the dimension of the space Sym(2n) of symmetric 
2n x 2 n matrices S = S T , given by dim Sym(2n) = ^(2n)(2n + 1) = n(2n+ 1). 
□ 

Remark 3 The space of Liouvillian forms with linear coefficients on a sym- 
plectic manifold (M, ui) is isomorphic to sp(M,u>) as an affine space. 

5.2 Liouvillian forms on the product manifold 

We are interested on the space of Liouvillian forms with linear components on 
the product manifold (M,coq) given by 

6 = nlOi - 7r*d 2 , 9 % er(T*tf i ),i= 1,2, ( 12 ) 

where 6 i,i = 1,2 are represented in local coordinates by 9 t = dz{\J + Si)z, 
i = 1,2 (see Lemma l54l . 

Lemma 5.4 Consider local coordinates (z,Z) £ (M,uq), where z G M\ and 
Z £ M 2 . Liouvillian forms on ( M,u)q ) given in 112 1) have a representation 
9 = (dz,dZ)A(z, Z) T for the matrix A = S + \ J , where A £ M 4 „x 4 n(R) is a 
symmetric matrix with the form 

S = ( ^ ) , Si, 5 2 ,0 2n G M 2 nx 2 n (R), S l = Sj, i = 1,2. 

Proof. Consider a point on the product manifold (z,Z) T £ ( M,uiq) written in 
local coordinates of the factors z £ M\ and Z £ M 2 . Using Lemma [5TT1 6 has a 
representation 


9 = ( dz,dZ)A 



A — S + 


(13) 


where S is a symmetric matrix of size 4n x 4 n. The key property is the decom¬ 
position of the tangent bundle TM = TM\ © TM 2 , since the representation of 
Liouvillian forms of type m is given by a matrix with two blocks 


( \J + Si 0 2n 

V “ 0 2n ~(\J + S 2 ) 


S 1 0 2n 
0 2n S‘2 


if J 0 2n 
2 'v 0 2n —J 


s + y 


proving the lemma. □ 

The (affine) space of Liouvillian forms with linear coefficients on (M, w©) has 
dimension dim Sym^n) = 2n(4n+ 1) and the dimension of the subset of forms 
given by m is 2 dim Sym(2n ) = 2n(2n+l). However, just a subset of this space 
is adapted for the construction of numerical schemes. We will be interested in 
the particular case S = Si = — S 2 whose projection satisfies Condition |T] We 
will return to this argument when we study Hamiltonian systems in Section [8] 
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6 The implicit map from the projection of the 
Lagrangian submanifold 

The linearization of 9 give us the local expresion of the vertical bundle^ UA C 
TM to the Lagrangian submanifold A and we want to recover the tangent 
bundle T Ac TM. Since the tangent bundle TM accepts a decomposition by 
TM = VA ®TA and the tangent bundle is mapped into the vertical bundle by 
14A = JT Z A for every z € A, then we recover locally the tangent spaces using 
the complex structure T z A = J T (V Z A). We project this subbundle with the 
tangent projection 

TN = Ttt(TA) = Ttt(J t (V A)). (14) 


Lemma 6.1 Let (M, Q\ x Q 2 , 6, 7 r, T) be a special symplectic manifold and A C 
M be an integral submanifold for 6 = tt* 9 i —tt %62 where 9i = dz(^J+Si)z. Then 
for every x £ A, the subspaces T ni ^Ai := Tni(J T (V X A)) C T x Mi, i = 1,2 are 
Lagrangian with local expression Li = (ll — JSi)z, satisfying J + Si)z = 0, 
i = 1,2. 

Proof. Applying point 2) from Lemma 12.11 the projection of the Lagrangian 
submanifold A £ M, by 7iy(A) = A* C Mj is a Lagrangian submanifold in Mi, 
i = 1,2. In local coordinates, the expressions 

Li = {\l — JS') z and (A J + Sf) z = 0, i = 1,2, (15) 

are the equations of the tangent space to the integral submanifold A* C Mi 
defined by the Liouvillian form 9i = dzQ^J + Sf)z. They define the same 
tangent space since A is an integral submanifold for 9 = n*9i — tt %92 and A j are 
integral submanifolds for 64 □ 

We define the linear space given by the sum V = T xl ( x )Ai + T„ 2 ^ A 2 as 
vector subspaces of TM. Using Lemma [6TT1 we can write V = L 1 +L 2 . We want 
that V satisfies Condition [l] however, point 2) in Condition |T| is enough for a 
symplectic map as we will see below. 

Lemma 6.2 Define the implicit map induced by the sum of the linear spaces 
V = L\ + L 2 from Lemma \6.1[ given by 

p(z, Z) = (\l 2 r 1 — JS 1 ) z + ( \l2n + JS 2 ) Z. (16) 

Then, p{z, z) = z, if and only if Si = S' 2 . 

Proof. Evaluating in (z, z) G A C M we have p(z, z) = z + J(5 < 2 — Si)z which 
produces the identity on the diagonal, if and only if J(S 2 — S±)z = 0. □ 

*The null space or the kernel of de projection Tit. 
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The implicit map induced by the projection of the lagrangian submanifold 
A adapted for the construction of symplectic maps is given in local coordinates 
of the product manifold (z, Z) G M by 

p{z,Z) = \{z + Z) + b{Z-z), (17) 

where b € M 2 nX 2 n(®.) is a Hamiltonian matrix b T J + Jb = 0. 


7 The method of Liouvillian forms 

The construction of implicit symplectic maps using the method of Liouvillian 
forms is obtained by restructuring the construction of the last section avoid¬ 
ing the explicit use of the special symplectic manifold. We recover all the in¬ 
formation from two different (but related) Liouvillian forms and the complex 
structure. 

1) For any exact symplectic manifold (M.uj) consider two copies (M t .oj 7 ), 
i = 1,2 for the construction of the product manifold (M,w@) and select a 
Liouvillian form 9 on M . 

We define the Liouvillian form satisfying Condition [I] in the following way: 
fix a Liouvillian form 9\ with linear coefficients on (Mi,wi) and express it in 
local coordinates by 9\ = dz(\J + S)z; define a second Liouvillian form on 
(M 2 ,u> 2 ) by 62 = dZ(^J — S)Z. The form 

9 = 7T*01 — 71202, (18) 

is Liouvillian on the product manifold. The expression of the vertical bundle in 
local coordinates is 

0 I s -ij"+s)(i)' se "“><«)■ s = sT 

2) The Pffafian equation 9 = 0 defines implicitly and locally the Lagrangian 
submanifold A. We use the complex structure J associated to the symplectic 
form uq of the product manifold to obtain the equations of the tangent spaces 
by T A = J T V A. In local coordinates we have 

t~" 0 ( z ) • 6 = JS e s = sT 

3) Project this Lagrangian subspace on the tangent bundle of the original 
symplectic manifold TM with Tni(TA) + Ttt 2 (T A) as linear subspaces. We will 
have the induced implicit map 

p{z,Z) = (\h n ~ b) z + (\h n + b) Z, (19) 

4) Verify if 




TM = Tiri(TA) ® Ttt 2 (TA) = T Ai ® TA 2 , 
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which holds when the projected Lagrangian submanifolds are complementary. 
Since z and Z are two different close points, we only check if the implicit map 
restricted to the diagonal is the identity p\a = l2n- The map m satisfies this 
condition by construction. 

Theorem 7.1 Let (M, w) be an exact symplectic manifold of dimension 2 n and 
U C M a start-shaped open subset containing the points z,Z G U. Define an 
implicit map p : U x U —> U by expression m • The implicit map (z, Z) K > 
Z = p(z, Z) is symplectic. 

Proof. We prove this result using implicit differentiation. Consider the implicit 
mapping 4' given by 


&(z,Z) = Z-p(z,Z) = 0. 


( 20 ) 


Implicit differentiation of (BUI) gives 


d&(z,Z) 

dz 


= -u. 


2 n 


+ b and 


Denoting the partial derivatives of T by 
&Sf(z,Z) 


B, = — 


dz 


d^{z,Z) 


dZ 


= bhn - b. 


( 21 ) 


and Bn = 


d^(z, Z) 
dZ : 


the amplification matrix of the linearized system mi is given by B = B 2 1 o B \, 
and T is symplectic if the matrix B is symplectic. We use the fact that the 
transpose B T = Bf o Bf T of a symplectic matrix is symplectic to tranform 
the symplecticity condition into Bf 1 o B\ JBf o Bf T = J, or equivalently into 
B\ JBi — B 2 JB 2 = 0. This condition is satisfied since B\ = B 2 , and the implicit 
map is symplectic as we want to prove. □ 


Remark 4 The converse is in general not true since we are characterizing only 
symplectic maps solving a variational problem with particular boundary condi¬ 
tions. 


We can consider the point z = p(z, Z) as an intermediate point on the path 
joining z to Z minimizing some variational problem. To develop this idea, we 
introduce some additional terminology. An implicit map <f> : U x U —> U is 
called consistent if there exist two explicit maps Vh > 4 > 2 '■ U —> U and a point 
z GU, such that 


z = tpi(z) and z = i/j 2 (Z). (22) 

We say that z is the point of consistency and if = ipf 1 o ip 1 : U U its 
consistency map. It is an explicit well defined map. We say that 4> interleaves 
a symplectic map if its consistency map is symplectic. 

We need a result from the Weyl’s extension of Cayley’s transformation. The 
interested readers are referred to m for the proof. 
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Lemma 7.2 (Generalized Cayley’s transformation) If the non-exceptional 
matrice ^ ] H and S are connected by the relations 

S = (/ — H)(7 + H) -1 = (I + H) _1 (7 — H) 

H = (/-S)(/ + S)- 1 = (7 + S)- 1 (/-S) 

and G is any matrix, then S T GS = G if and only ifH T G + GH = 0. 

For the symplectic case, we fix the arbitrary matrix G = J to the complex 
structure on TM. Now we can relate the induced implicit map p with Cayley’s 
transformation in the following 

Proposition 2 With the same hypotheses as Theorem \7.1\ consider the explicit 
map ip = ip 2 1 oip 1 : U ~t U as the consistency map associated to the implicit map 
p : TU x TU —>• TU given in 119\) . Then, the consistency map ip is symplectic 
and it corresponds to Cayley’s transformation of the matrix 2 JS. 

Proof. Consider p(z, Z) as a linear combination of two explicit linear maps com¬ 
ing from z and Z, in the form p(z, Z) = ^ (Tipi( z ) + Tip 2 (Z)). From expression 
m , we can write explicitely 

T1> 1 (z) = (I 2n -2JS)z and Tip 2 {Z) = {I 2n + 2JS)Z. (23) 

Since S E M 2nX 2 n(f&) is a symmetric matrix, the matrix H = 2 JS is Hamil¬ 
tonian. In this case, the consistency map associated to p(z, Z) is given by 
ip = tp 2 1 o ipi, whose linearization gives 

Tip = TiPi 1 o Tipi = (hn + H )-\h n - H). (24) 

This is the symplectic Cayley’s transformation of the Hamiltonian matrix H. 
Applying Weyl’s Lemma we obtain immediately that the consistency map ip : 
U B z Z = ip(z) € U is symplectic. □ 

Corollary 3 For every symplectic mapping ip E Sp(M, w) on an exact sym¬ 
plectic manifold, there exists a Liouvillian form 9^ adapted to ip. 

Proof. Suppose ip E Sp(M , w) is represented at every point m E M by a non- 
exceptional matrix S E Sp(2n). There exists an associated Hamiltonian matrix 
given by the symplectic Cayley’s transformation H = (I 2n — S )(/ 2n + S) -1 . 
Consequently, the linearization of the Liouvillian form 9^ adapted to ip is given 
by some square matrix R = i J ( I 2n + H) in the form 9^ = dzRz. Moreover, 
the implicit map p : TU x TU —> TU in coordinates (z,Z) E U x U is defined 
in terms of H by 

P(z,Z) = (ll 2n -lHjz+(^I 2n + ^Hjz. (25) 

6 A matrix A E GL(n) is said to be non-exceptional if det (I + A) / 0, where I is the 
identity matrix in GL(n). 
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□ 

Consider the matrix G in Weyl’s lemma [7721 as a bilinear form u T Gv. The 
trivial case G = l 2 n gives the Euclidean metric and the Cayley’s transformation 
relates symmetric with antisymmetric matrices. The case G = J is the symplec- 
tic form and it relates symplectic with Hamiltonian matrices. In JBj we show 
that this property with an additional dependency on a parameter discretizing 
time produce reversible maps. 


8 Hamiltonian systems and implicit symplectic 
integrators 

The numerical simulation of Hamiltonian vector fields is made by numerical 
schemes called symplectic algorithms or symplectic integrators. We introduce the 
notation and terminology for the application of Liouvillian forms in symplectic 
integration. 

A Hamiltonian system (M, ui, Xh) is a vector field X = Xh on a symplectic 
manifold (M, ui) such that inner product of the symplectic form and the vector 
field satisfies ix H w = ~dH, for a differentiable function H : M —> R. 

The set of solutions : M x R —>• M of the Hamiltonian vector field Xh 
is called the flow of Xh and it is defined by 

jy H {z) = X H {p t H {z)). (26) 

The flow is a 1-parameter subgroup of symplectic diffeomorphisms with param¬ 
eter t £ R. For every point z 0 £ M, the solution Ph(zo) is the integral curve of 
Xh passing by zo at time t = 0. 

A symplectic algorithm with stepsize h, is the numerical approximation ifh 
of the map for t = h fixed, which is smooth with respect to h and H , and 
preserves the symplectic form (? fh)*w = w. Consider an open set U £ M and 
two points on the flow of Xh, say Zt = <p*{zo) and Zt+h = <PH +k \ z o)- 
group property of the flow, it is enough to perform the analysis around t = 0, 
then the points will be denoted by zq and Zh- 

Theorem 8.1 Let ( M,9,Xh ) be a Hamiltonian system on an exact symplectic 
manifold. Consider a convex open set U C M containing the points z o and Zh 
on the flow of Xh ■ If an implicit map p : U x U —* U is defined by 

p(z 0 ,z h ) = (\h n - b) z 0 + (\h n + b) z h , (27) 

where b is a Hamiltonian matrix in Af 2 ra x 2 n (®1) ; then, the map 

Zh = z 0 + hX H o p(z 0 , z h ) 


is symplectic. 
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Proof. Since Xh is a Hamiltonian vector field, invariant under symplectic trans¬ 
formations, it suffices that p(zo,Zh) be a symplectic map. Applying Theorem 
o we obtain the desired result. □ 

An alternative proof of the theorem is using implicit differentiation. Consider 
the implicit mapping T given by 


^(z 0 ,z h ) = z h - z 0 - hX H o p(z 0 ,z h ) = 0 


(28) 


as in the proof of Teorem 17.11 Implicit differentiation of (E51) using the chain 
rule gives 


B 1 :=- 
B 2 := 


d'S'(z 0 ,z h ) 

dz 0 

dV(zo,z h ) 

dz h 


hn — hH ( \l 2n — b) 
—hn — hH (\l 2n + b) 


(29) 

(30) 


where H is a Hamiltonian matrix given by H = JD 2 H and D 2 H is the Hessian 
matrix of H. The symplecticity test BiJBf — B 2 JBj = 0 is satisfied since H 
and b are Hamiltonian matrices. 


8.1 The generalized implicit symplectic Euler scheme 

The relevance of this method is that we have a linear (continuous) space of di¬ 
mension n(2n+l) where we can select Hamiltonian matrices for the construction 
of numerical schemes. Using Theorem 18.II we are able to define the generalized 
implicit symplectic Euler scheme as the map given by 

iph ■ U x U H- U 

(z 0 , z h ) i->- z h = z 0 + hX H ( z) 

where z = ~(zo + Zh) + b{zh — ~o) and it corresponds to the map z 
given by the expression m- 

Remark 5 In a slightly different context, Feng Kang (unpublished) has shown 
that generating functions obtained by some particular type of matrices a £ 
Af 4 „ X 4 rt(R) can be reduced to generating functions constructed by an equiva¬ 
lent simplified matrix containing the two submatrices given in expression m 
(see Ge and Dau-liu In their method, those matrices are the input for the 
evolutive Hamilton-Jacobi equation. 

Proposition 3 The elementary Liouvillian form on an exact symplectic 
manifold (M,u>), produces a Liouvillian form on the product manifold ( M,u>q ) 
given by 6 = tt(9o — n^Oo = IJzqJzq — IdzhJzh- In Darboux’s coordinates it 
corresponds to 


(31) 
= p(zo,Zh) 


9 = \ (podqo ~ qodpo - Phdqh + qhdqn ), 


(32) 


whose generalized symplectic Euler scheme is the symplectic mid-point rule. 
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Proof. It is a direct computation using the fact that elementary Liouvillian form 
does not have a symmetric component. In terms of the generalized symplectic 
Euler scheme, their argument z have a null Hamiltonian matrix b = 02 ra . The 
argument z = i(z 0 + Zh) gives the scheme 

Zh = z 0 + hX H (i(z 0 + z h )) , (33) 

which is the mid-point rule. □ 

Proposition 4 In terms of the generalized symplectic Euler scheme m with 
Darboux’s coordinates zo = ( qo,Po) T and Zh = (qh,Ph) T , the Hamiltonian ma¬ 
trices 

&/1 = 5 ( 5), = b G M 2 nx2n(R). 

corresponds to the symplectic Euler A and B schemes, respectively. 

Proof. Write the sum and difference vectors in Darboux’s coordinates by 

Computing b A (z h - zo) leads to b A (z h - z 0 ) = \ (qo - qh,Ph ~ Po) T and conse¬ 
quently z = i(zo + Zh) + b A (zh — zq) produces z = {qo,Ph) T ■ In the same way 
z = i(z 0 + Zh) + bsizh — zo) produces z = (qh,Po) T - Using these arguments in 
the generalized scheme m we obtain the standard Euler schemes. □ 

9 Two families of implicit symplectic integrators 

In this section we develop two examples of the application of the method of 
Liouvillian forms for constructing symplectic integrators. The first one consid¬ 
ers the construction of the special symplectic manifold on Q\ x Q 2 , using the 
symplectomorphism T : M —> T*(Qi x <2 2 )- The second one is given for apply¬ 
ing the method directly. In the practice, we will use the generalized symplectic 
Euler scheme given in (1311) . 

9.1 Example 1: A simple family from symplectic rotations 

Here, we explicitely define the symplectomorphism between the product mani¬ 
fold (M,u>e) and the cotangent bundle T*(Qi x Q 2 ) using a family of symplec¬ 
tic rotations. We take the loop of symplectic rotations induced by the complex 
structure J i —> exp (<j)J) on the product manifold (M,ljq). 

Lemma 9.1 If J is a complex structure then exp (<j>J) = cos((/>)/ 4 n + sin {(f) J. 
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Proof. Use J 2 = — R n in the development of the exponential. □ Define the 
symplectic rotation R 0 € Sp(M,ojq) by = exp (</>J) and define the curve of 
diffeomorphisms ^ = E\ o i? 0 , cj) g [ 0 , 7 r/ 2 ]. 

Given symplectic coordinates {qo,Po, Qh,Ph) G M and (x 0 , x h , y 0 , yh) G T*(Qi x 
Q 2 ), the diffeomorphism 4/0 : M —> T*(Qi x <2 2 ) pulls-back the Liouville form 
9q 1 xQ 2 = Vodxo + yhdxh into 9 0 = 4/^ 9q iX q 2 expressed in these coordinates 
by 

00 = cos 2 (fpodqo - sin 2 (fph.dqh - sin 2 4>q 0 dp 0 + cos 2 (fqndph 

+ cos <f> sin <p(qodqo - p 0 dp 0 ) + cos </>sin (j)(q h dqh - Ph.dph) (34) 

The quintuple (Qi x Q 2 ), 00 , 7 r, 4 / 0 j is a family of special symplectic man¬ 
ifolds on Qi x Q 2 . We have the following 

Lemma 9.2 The set of forms \34\) is a 1-parameter family of Liouvillian forms. 
The path 00 contains the Liouvillian forms associated to the generating functions 
of type II and type III. 

Proof. By a direct computation we obtain g ?00 = dpo A dqo — dph A dqh = oJq, 
showing that any element in the family is Liouvillian on The forms at 

the values (f = 0 and <j> = j are 

00 ( 0=0 = Podqo + qhdph and 00 | 0 =7r /2 = -qodpo ~ Phdqh- 

By Lemma 14.21 they are associated to generating functions of type II and III , 
respectively. □ 

The family of vertical fibers to A 0 C M, associated to the Liouvillian forms 
00 is given by equations 

< 7 o = cos 2 tppo + cos <f> sin (f>q 0 , Po = ~ sin 2 0g o — cos 4>i sin <pp 0 

Qh = - sin 2 (fph + cos <j>i sin 4>q h , p h = cos 2 f>q h - cos <j) sin 4>p h . 

We obtain the tangent fibers by left multiplying vertical fibers by J T , map¬ 
ping (qo,Po,q h ,Ph) i->- {-po,qo,Ph,~qh)- Finally, the projection Tir^A^) ^ 
T 7 ri(TA 0 )®T 7 r 2 (TA 0 ) in local coordinates is p(q 0 ,p 0 , qh,Ph) = {Ph~Po,qo ~ qh), 
or explicitely by 

qh = cos 2 <j>q h + sin 2 4 >q 0 + cos f> sin <j)(p 0 - ph) 

Ph = cos 2 4>p 0 + sin 2 (j)p h + cos sin 4>(q 0 - q h ). (35) 

The family of implicit symplectic integrators is given by 

9h = qo + | (■ q h ,p h ) Ph=Po~ | ( q h ,Ph ) • (36) 

Evaluating these coordinates in <fi = 0 and <f> = 7t/ 2 we obtain 

{Qh,Ph)\^ 0 = (Qh,Po) and {qh,Ph)\ (j> = v /2 = (Qo,Ph) • (37) 

We have proved the following result. 


23 


Corollary 4 The implicit scheme Zh = zo+h.X H op(zQ, Zh) where the linear map 
p(zo,Zh) = (q h ,ph) is defined by the expression H35\ ) , is a family of symplectic 
integrators which joins the symplectic Euler schemes A and B. 


Using the identities 


cos <f> sin (f = } sin 2 cf> = cos 2 0 = 1+c ° s2 ^ , 

we rewrite (l35l) in compact form z = (qh,Ph) by z = \{zq + Zh) + b{zh — zq), 
where 


b = 


1 

2 


— cos 2<f> sin 2<f> \ 
sin 2 (f> cos 2</> J ’ 


b £ M2,T,x2n(R)- 


(38) 


9.2 Example 2: A three-parameter family 

The method of Liouvillian forms applied to the construction of symplectic inte¬ 
grators yields the generalized symplectic Euler scheme (1311) . This scheme uses 
the argument z = \{zq + Zh) + b(zh — zq), where b £ M 2 X 2(®0 is an arbitrary 
Hamiltonian matrix. In a problem with 1 degree of freedom, b is a matrix of 
the form 

*=(“ _? a ). a,^£R. 

These are three free parameters in the generalized Euler method. Explicitely 
we have the numerical scheme 

qh=qo + h^-(q,p) and p h = Po ~ (q,p), (39) 

where 

q = (I - of) q Q + (i + a) q h + /3{p h - p 0 ) 

P = (5 + a) q 0 + (5 - a) q h + j(qh - qo), 

and a, fS, 7 < h 2 < ^ have small values. Parameter a modifies the symmetry of 
the oscillations with respect to time and the limiting cases (a, /3, 7 ) = (— 0,0) 

and (a,/3, 7 ) = (^,0,0) are the usual symplectic Euler A and B integrators. 

The parameters {a, /3, 7 } modify the error and the oscillations of the nu¬ 
merical integrator around the orbits of the Hamiltonian flow. This fact has a 
geometrical justification that we develop in [5]. In a parallel article [ 6 ], we make 
a numerical study of this family explaining the oscillations of the numerical 
solution around the exact solution using a variational point of view. 
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